Dysregulation of CD4+ and CD8+ resident memory T, myeloid, and stromal cells in steroid-experienced, checkpoint inhibitor colitis.

BACKGROUND: Colitis caused by checkpoint inhibitors (CPI) is frequent and is treated with empiric steroids, but CPI colitis mechanisms in steroid-experienced or refractory disease are unclear. METHODS: Using colon biopsies and blood from predominantly steroid-experienced CPI colitis patients, we performed multiplexed single-cell transcriptomics and proteomics to nominate contributing populations. RESULTS: CPI colitis biopsies showed enrichment of CD4+resident memory (RM) T cells in addition to CD8+ RM and cytotoxic CD8+ T cells. Matching T cell receptor (TCR) clonotypes suggested that both RMs are progenitors that yield cytotoxic effectors. Activated, CD38+ HLA-DR+ CD4+ RM and cytotoxic CD8+ T cells were enriched in steroid-experienced and a validation data set of steroid-naïve CPI colitis, underscoring their pathogenic potential across steroid exposure. Distinct from ulcerative colitis, CPI colitis exhibited perturbed stromal metabolism (NAD+, tryptophan) impacting epithelial survival and inflammation. Endothelial cells in CPI colitis after anti-TNF and anti-cytotoxic T-lymphocyte-associated antigen 4 (anti-CTLA-4) upregulated the integrin α4ß7 ligand molecular vascular addressin cell adhesion molecule 1 (MAdCAM-1), which may preferentially respond to vedolizumab (anti-α4ß7). CONCLUSIONS: These findings nominate CD4+ RM and MAdCAM-1+ endothelial cells for targeting in specific subsets of CPI colitis patients.

Defining and benchmarking open problems in single-cell analysis.

With the growing number of single-cell analysis tools, benchmarks are increasingly important to guide analysis and method development. However, a lack of standardisation and extensibility in current benchmarks limits their usability, longevity, and relevance to the community. We present Open Problems, a living, extensible, community-guided benchmarking platform including 10 current single-cell tasks that we envision will raise standards for the selection, evaluation, and development of methods in single-cell analysis.

Single-cell and spatial multi-omics highlight effects of anti-integrin therapy across cellular compartments in ulcerative colitis.

Ulcerative colitis (UC) is driven by immune and stromal subsets, culminating in epithelial injury. Vedolizumab (VDZ) is an anti-integrin antibody that is effective for treating UC. VDZ is known to inhibit lymphocyte trafficking to the intestine, but its broader effects on other cell subsets are less defined. To identify the inflammatory cells that contribute to colitis and are affected by VDZ, we perform single-cell transcriptomic and proteomic analyses of peripheral blood and colonic biopsies in healthy controls and patients with UC on VDZ or other therapies. Here we show that VDZ treatment is associated with alterations in circulating and tissue mononuclear phagocyte (MNP) subsets, along with modest shifts in lymphocytes. Spatial multi-omics of formalin-fixed biopsies demonstrates trends towards increased abundance and proximity of MNP and fibroblast subsets in active colitis. Spatial transcriptomics of archived specimens pre-treatment identifies epithelial-, MNP-, and fibroblast-enriched genes related to VDZ responsiveness, highlighting important roles for these subsets in UC.

Single-cell and spatial multi-omics highlight effects of anti-integrin therapy across cellular compartments in ulcerative colitis.

Ulcerative colitis (UC) is driven by immune and stromal subsets, culminating in epithelial injury. Vedolizumab (VDZ) is an anti-integrin antibody that is effective for treating UC. VDZ is known to inhibit lymphocyte trafficking to the intestine, but its broader effects on other cell subsets are less defined. To identify the inflammatory cells that contribute to colitis and are affected by VDZ, we performed single-cell transcriptomic and proteomic analyses of peripheral blood and colonic biopsies in healthy controls and patients with UC on VDZ or other therapies. Here we show that VDZ treatment is associated with alterations in circulating and tissue mononuclear phagocyte (MNP) subsets, along with modest shifts in lymphocytes. Spatial multi-omics of formalin-fixed biopsies demonstrates trends towards increased abundance and proximity of MNP and fibroblast subsets in active colitis. Spatial transcriptomics of archived specimens pre-treatment identifies epithelial-, MNP-, and fibroblast-enriched genes related to VDZ responsiveness, highlighting important roles for these subsets in UC.

scDiffCom: a tool for differential analysis of cell-cell interactions provides a mouse atlas of aging changes in intercellular communication.

Dysregulation of intercellular communication is a hallmark of aging. To better quantify and explore changes in intercellular communication, we present scDiffCom and scAgeCom. scDiffCom is an R package, relying on approximately 5,000 curated ligand-receptor interactions, that performs differential intercellular communication analysis between two conditions from single-cell transcriptomics data. Built upon scDiffCom, scAgeCom is an atlas of age-related cell-cell communication changes covering 23 mouse tissues from 58 single-cell RNA sequencing datasets from Tabula Muris Senis and the Calico murine aging cell atlas. It offers a comprehensive resource of tissue-specific and sex-specific aging dysregulations and highlights age-related intercellular communication changes widespread across the whole body, such as the upregulation of immune system processes and inflammation, the downregulation of developmental processes, angiogenesis and extracellular matrix organization and the deregulation of lipid metabolism. Our analysis emphasizes the relevance of the specific ligands, receptors and cell types regulating these processes. The atlas is available online ( https://scagecom.org ).

Developmental origin and local signals cooperate to determine septal astrocyte identity.

Astrocyte specification during development is influenced by both intrinsic and extrinsic factors, but the precise contribution of each remains poorly understood. Here we show that septal astrocytes from Nkx2.1 and Zic4 expressing progenitor zones are allocated into non-overlapping domains of the medial (MS) and lateral septal nuclei (LS) respectively. Astrocytes in these areas exhibit distinctive molecular and morphological features tailored to the unique cellular and synaptic circuit environment of each nucleus. Using single-nucleus (sn) RNA sequencing, we trace the developmental trajectories of cells in the septum and find that neurons and astrocytes undergo region and developmental stage-specific local cell-cell interactions. We show that expression of the classic morphogens Sonic hedgehog (Shh) and Fibroblast growth factors (Fgfs) by MS and LS neurons respectively, functions to promote the molecular specification of local astrocytes in each region. Finally, using heterotopic cell transplantation, we show that both morphological and molecular specifications of septal astrocytes are highly dependent on the local microenvironment, regardless of developmental origins. Our data highlights the complex interplay between intrinsic and extrinsic factors shaping astrocyte identities and illustrates the importance of the local environment in determining astrocyte functional specialization.

Single-nuclei characterization of pervasive transcriptional signatures across organs in response to COVID-19.

Background: Infection by coronavirus SARS-CoV2 is a severe and often deadly disease that has implications for the respiratory system and multiple organs across the human body. While the effects in the lung have been extensively studied, less is known about the impact COVID-19 has across other organs. Methods: Here, we contribute a single-nuclei RNA-sequencing atlas comprising six human organs across 20 autopsies where we analyzed the transcriptional changes due to COVID-19 in multiple cell types. The integration of data from multiple organs enabled the identification of systemic transcriptional changes. Results: Computational cross-organ analysis for endothelial cells and macrophages identified systemic transcriptional changes in these cell types in COVID-19 samples. In addition, analysis of gene modules showed enrichment of specific signaling pathways across multiple organs in COVID-19 autopsies. Conclusions: Altogether, the COVID Tissue Atlas enables the investigation of both cell type-specific and cross-organ transcriptional responses to COVID-19, providing insights into the molecular networks affected by the disease and highlighting novel potential targets for therapies and drug development. Funding: The Chan-Zuckerberg Initiative, The Chan-Zuckerberg Biohub.

Aging Fly Cell Atlas identifies exhaustive aging features at cellular resolution.

Aging is characterized by a decline in tissue function, but the underlying changes at cellular resolution across the organism remain unclear. Here, we present the Aging Fly Cell Atlas, a single-nucleus transcriptomic map of the whole aging Drosophila. We characterized 163 distinct cell types and performed an in-depth analysis of changes in tissue cell composition, gene expression, and cell identities. We further developed aging clock models to predict fly age and show that ribosomal gene expression is a conserved predictive factor for age. Combining all aging features, we find distinctive cell type-specific aging patterns. This atlas provides a valuable resource for studying fundamental principles of aging in complex organisms.

The mitochondrial unfolded protein response regulates hippocampal neural stem cell aging.

Aging results in a decline in neural stem cells (NSCs), neurogenesis, and cognitive function, and evidence is emerging to demonstrate disrupted adult neurogenesis in the hippocampus of patients with several neurodegenerative disorders. Here, single-cell RNA sequencing of the dentate gyrus of young and old mice shows that the mitochondrial protein folding stress is prominent in activated NSCs/neural progenitors (NPCs) among the neurogenic niche, and it increases with aging accompanying dysregulated cell cycle and mitochondrial activity in activated NSCs/NPCs in the dentate gyrus. Increasing mitochondrial protein folding stress results in compromised NSC maintenance and reduced neurogenesis in the dentate gyrus, neural hyperactivity, and impaired cognitive function. Reducing mitochondrial protein folding stress in the dentate gyrus of old mice improves neurogenesis and cognitive function. These results establish the mitochondrial protein folding stress as a driver of NSC aging and suggest approaches to improve aging-associated cognitive decline.

Indian Hedgehog release from TNF-activated renal epithelia drives local and remote organ fibrosis.

Progressive fibrosis is a feature of aging and chronic tissue injury in multiple organs, including the kidney and heart. Glioma-associated oncogene 1 expressing (Gli1+) cells are a major source of activated fibroblasts in multiple organs, but the links between injury, inflammation, and Gli1+ cell expansion and tissue fibrosis remain incompletely understood. We demonstrated that leukocyte-derived tumor necrosis factor (TNF) promoted Gli1+ cell proliferation and cardiorenal fibrosis through induction and release of Indian Hedgehog (IHH) from renal epithelial cells. Using single-cell-resolution transcriptomic analysis, we identified an "inflammatory" proximal tubular epithelial (iPT) population contributing to TNF- and nuclear factor κB (NF-κB)-induced IHH production in vivo. TNF-induced Ubiquitin D (Ubd) expression was observed in human proximal tubular cells in vitro and during murine and human renal disease and aging. Studies using pharmacological and conditional genetic ablation of TNF-induced IHH signaling revealed that IHH activated canonical Hedgehog signaling in Gli1+ cells, which led to their activation, proliferation, and fibrosis within the injured and aging kidney and heart. These changes were inhibited in mice by Ihh deletion in Pax8-expressing cells or by pharmacological blockade of TNF, NF-κB, or Gli1 signaling. Increased amounts of circulating IHH were associated with loss of renal function and higher rates of cardiovascular disease in patients with chronic kidney disease. Thus, IHH connects leukocyte activation to Gli1+ cell expansion and represents a potential target for therapies to inhibit inflammation-induced fibrosis.

Concordance of MERFISH spatial transcriptomics with bulk and single-cell RNA sequencing.

Spatial transcriptomics extends single-cell RNA sequencing (scRNA-seq) by providing spatial context for cell type identification and analysis. Imaging-based spatial technologies such as multiplexed error-robust fluorescence in situ hybridization (MERFISH) can achieve single-cell resolution, directly mapping single-cell identities to spatial positions. MERFISH produces a different data type than scRNA-seq, and a technical comparison between the two modalities is necessary to ascertain how to best integrate them. We performed MERFISH on the mouse liver and kidney and compared the resulting bulk and single-cell RNA statistics with those from the Tabula Muris Senis cell atlas and from two Visium datasets. MERFISH quantitatively reproduced the bulk RNA-seq and scRNA-seq results with improvements in overall dropout rates and sensitivity. Finally, we found that MERFISH independently resolved distinct cell types and spatial structure in both the liver and kidney. Computational integration with the Tabula Muris Senis atlas did not enhance these results. We conclude that MERFISH provides a quantitatively comparable method for single-cell gene expression and can identify cell types without the need for computational integration with scRNA-seq atlases.

A 2-Gene Host Signature for Improved Accuracy of COVID-19 Diagnosis Agnostic to Viral Variants.

The continued emergence of SARS-CoV-2 variants is one of several factors that may cause false-negative viral PCR test results. Such tests are also susceptible to false-positive results due to trace contamination from high viral titer samples. Host immune response markers provide an orthogonal indication of infection that can mitigate these concerns when combined with direct viral detection. Here, we leverage nasopharyngeal swab RNA-seq data from patients with COVID-19, other viral acute respiratory illnesses, and nonviral conditions (n = 318) to develop support vector machine classifiers that rely on a parsimonious 2-gene host signature to diagnose COVID-19. We find that optimal classifiers include an interferon-stimulated gene that is strongly induced in COVID-19 compared with nonviral conditions, such as IFI6, and a second immune-response gene that is more strongly induced in other viral infections, such as GBP5. The IFI6+GBP5 classifier achieves an area under the receiver operating characteristic curve (AUC) greater than 0.9 when evaluated on an independent RNA-seq cohort (n = 553). We further provide proof-of-concept demonstration that the classifier can be implemented in a clinically relevant RT-qPCR assay. Finally, we show that its performance is robust across common SARS-CoV-2 variants and is unaffected by cross-contamination, demonstrating its utility for improved accuracy of COVID-19 diagnostics. IMPORTANCE In this work, we study upper respiratory tract gene expression to develop and validate a 2-gene host-based COVID-19 diagnostic classifier and then demonstrate its implementation in a clinically practical qPCR assay. We find that the host classifier has utility for mitigating false-negative results, for example due to SARS-CoV-2 variants harboring mutations at primer target sites, and for mitigating false-positive viral PCR results due to laboratory cross-contamination. Both types of error carry serious consequences of either unrecognized viral transmission or unnecessary isolation and contact tracing. This work is directly relevant to the ongoing COVID-19 pandemic given the continued emergence of viral variants and the continued challenges of false-positive PCR assays. It also suggests the feasibility of pan-respiratory virus host-based diagnostics that would have value in congregate settings, such as hospitals and nursing homes, where unrecognized respiratory viral transmission is of particular concern.

The hepatic integrated stress response suppresses the somatotroph axis to control liver damage in nonalcoholic fatty liver disease.

Nonalcoholic fatty liver disease (NAFLD) can be ameliorated by calorie restriction, which leads to the suppressed somatotroph axis. Paradoxically, the suppressed somatotroph axis is associated with patients with NAFLD and is correlated with the severity of fibrosis. How the somatotroph axis becomes dysregulated and whether the repressed somatotroph axis impacts liver damage during the progression of NAFLD are unclear. Here, we identify a regulatory branch of the hepatic integrated stress response (ISR), which represses the somatotroph axis in hepatocytes through ATF3, resulting in enhanced cell survival and reduced cell proliferation. In mouse models of NAFLD, the ISR represses the somatotroph axis, leading to reduced apoptosis and inflammation but decreased hepatocyte proliferation and exacerbated fibrosis in the liver. NAD+ repletion reduces the ISR, rescues the dysregulated somatotroph axis, and alleviates NAFLD. These results establish that the hepatic ISR suppresses the somatotroph axis to control cell fate decisions and liver damage in NAFLD.

Integrated host-microbe plasma metagenomics for sepsis diagnosis in a prospective cohort of critically ill adults.

We carried out integrated host and pathogen metagenomic RNA and DNA next generation sequencing (mNGS) of whole blood (n = 221) and plasma (n = 138) from critically ill patients following hospital admission. We assigned patients into sepsis groups on the basis of clinical and microbiological criteria. From whole-blood gene expression data, we distinguished patients with sepsis from patients with non-infectious systemic inflammatory conditions using a trained bagged support vector machine (bSVM) classifier (area under the receiver operating characteristic curve (AUC) = 0.81 in the training set; AUC = 0.82 in a held-out validation set). Plasma RNA also yielded a transcriptional signature of sepsis with several genes previously reported as sepsis biomarkers, and a bSVM sepsis diagnostic classifier (AUC = 0.97 training set; AUC = 0.77 validation set). Pathogen detection performance of plasma mNGS varied on the basis of pathogen and site of infection. To improve detection of virus, we developed a secondary transcriptomic classifier (AUC = 0.94 training set; AUC = 0.96 validation set). We combined host and microbial features to develop an integrated sepsis diagnostic model that identified 99% of microbiologically confirmed sepsis cases, and predicted sepsis in 74% of suspected and 89% of indeterminate sepsis cases. In summary, we suggest that integrating host transcriptional profiling and broad-range metagenomic pathogen detection from nucleic acid is a promising tool for sepsis diagnosis.

Polygenic enrichment distinguishes disease associations of individual cells in single-cell RNA-seq data.

Single-cell RNA sequencing (scRNA-seq) provides unique insights into the pathology and cellular origin of disease. We introduce single-cell disease relevance score (scDRS), an approach that links scRNA-seq with polygenic disease risk at single-cell resolution, independent of annotated cell types. scDRS identifies cells exhibiting excess expression across disease-associated genes implicated by genome-wide association studies (GWASs). We applied scDRS to 74 diseases/traits and 1.3 million single-cell gene-expression profiles across 31 tissues/organs. Cell-type-level results broadly recapitulated known cell-type-disease associations. Individual-cell-level results identified subpopulations of disease-associated cells not captured by existing cell-type labels, including T cell subpopulations associated with inflammatory bowel disease, partially characterized by their effector-like states; neuron subpopulations associated with schizophrenia, partially characterized by their spatial locations; and hepatocyte subpopulations associated with triglyceride levels, partially characterized by their higher ploidy levels. Genes whose expression was correlated with the scDRS score across cells (reflecting coexpression with GWAS disease-associated genes) were strongly enriched for gold-standard drug target and Mendelian disease genes.

Upper airway gene expression shows a more robust adaptive immune response to SARS-CoV-2 in children.

Unlike other respiratory viruses, SARS-CoV-2 disproportionately causes severe disease in older adults whereas disease burden in children is lower. To investigate whether differences in the upper airway immune response may contribute to this disparity, we compare nasopharyngeal gene expression in 83 children (<19-years-old; 38 with SARS-CoV-2, 11 with other respiratory viruses, 34 with no virus) and 154 older adults (>40-years-old; 45 with SARS-CoV-2, 28 with other respiratory viruses, 81 with no virus). Expression of interferon-stimulated genes is robustly activated in both children and adults with SARS-CoV-2 infection compared to the respective non-viral groups, with only subtle distinctions. Children, however, demonstrate markedly greater upregulation of pathways related to B cell and T cell activation and proinflammatory cytokine signaling, including response to TNF and production of IFNγ, IL-2 and IL-4. Cell type deconvolution confirms greater recruitment of B cells, and to a lesser degree macrophages, to the upper airway of children. Only children exhibit a decrease in proportions of ciliated cells, among the primary targets of SARS-CoV-2, upon infection. These findings demonstrate that children elicit a more robust innate and especially adaptive immune response to SARS-CoV-2 in the upper airway that likely contributes to their protection from severe disease in the lower airway.

The Tabula Sapiens: A multiple-organ, single-cell transcriptomic atlas of humans.

Molecular characterization of cell types using single-cell transcriptome sequencing is revolutionizing cell biology and enabling new insights into the physiology of human organs. We created a human reference atlas comprising nearly 500,000 cells from 24 different tissues and organs, many from the same donor. This atlas enabled molecular characterization of more than 400 cell types, their distribution across tissues, and tissue-specific variation in gene expression. Using multiple tissues from a single donor enabled identification of the clonal distribution of T cells between tissues, identification of the tissue-specific mutation rate in B cells, and analysis of the cell cycle state and proliferative potential of shared cell types across tissues. Cell type-specific RNA splicing was discovered and analyzed across tissues within an individual.

Dynamical Systems Model of RNA Velocity Improves Inference of Single-cell Trajectory, Pseudo-time and Gene Regulation.

Recent development in inferring RNA velocity from single-cell RNA-seq opens up exciting new vista into developmental lineage and cellular dynamics. However, the estimated velocity only gives a snapshot of how the transcriptome instantaneously changes in individual cells, and it does not provide quantitative predictions and insights about the whole system. In this work, we develop RNA-ODE, a principled computational framework that extends RNA velocity to quantify systems level dynamics and improve single-cell data analysis. We model the gene expression dynamics by an ordinary differential equation (ODE) based formalism. Given a snapshot of gene expression at one time, RNA-ODE is able to predict and extrapolate the expression trajectory of each cell by solving the dynamic equations. Systematic experiments on simulations and on new data from developing brain demonstrate that RNA-ODE substantially improves many aspects of standard single-cell analysis. By leveraging temporal dynamics, RNA-ODE more accurately estimates cell state lineage and pseudo-time compared to previous state-of-the-art methods. It also infers gene regulatory networks and identifies influential genes whose expression changes can decide cell fate. We expect RNA-ODE to be a Swiss army knife that aids many facets of single-cell RNA-seq analysis.

Transcriptomic profiling of blood from autoimmune hepatitis patients reveals potential mechanisms with implications for management.

Autoimmune hepatitis (AIH) is a poorly understood, chronic disease, for which corticosteroids are still the mainstay of therapy and most patients undergo liver biopsy to obtain a diagnosis. We aimed to determine if there was a transcriptomic signature of AIH in the peripheral blood and investigate underlying biologic pathways revealed by gene expression analysis. Whole blood RNA from 75 AIH patients and 25 healthy volunteers was extracted and sequenced. Differential gene expression analysis revealed 249 genes that were significantly differentially expressed in AIH patients compared to controls. Using a random forest algorithm, we determined that less than 10 genes were sufficient to differentiate the two groups in our cohort. Interferon signaling was more active in AIH samples compared to controls, regardless of treatment status. Pegivirus sequences were detected in five AIH samples and 1 healthy sample. The gene expression data and clinical metadata were used to determine 12 genes that were significantly associated with advanced fibrosis in AIH. AIH patients with a partial response to therapy demonstrated decreased evidence of a CD8+ T cell gene expression signal. These findings represent progress in understanding a disease in need of better tests, therapies, and biomarkers.